WAVELETS ON IRREGULAR GRIDS WITH ARBITRARY 
DILATION MATRICES, AND FRAME ATOMS FOR L 2 (R d ) 
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Abstract. In this article, wc develop a general method for constructing wa- 
velets {|det A j \ 1 / 2 7p(A j x — x j • j. ) : j £ J, k £ K} on irregular lattices of the 
form X = {xj t k £ R d : j £ J,k £ A"}, and with an arbitrary countable 
family of invertible d X d matrices {Aj £ GLd(K) : j £ J} that do not nec- 
essarily have a group structure. This wavelet construction is a particular case 
of general atomic frame decompositions of L 2 (R d ) developed in this article, 
that allow other time frequency decompositions such as non-harmonic Gabor 
frames with non-uniform covering of the Euclidean space lR d . Possible applica- 
tions include image and video compression, speech coding, image and digital 
data transmission, image analysis, estimations and detection, and seismology. 



1. Introduction 

Recently there has been a considerable interest in trying to obtain atomic decom- 
positions of the space L 2 (R d ). These decompositions are usually obtained in terms 
of frames generated by a family of functions translated on a regular grid, and di- 
lated by powers of a dilation matrix. The uniformity of the grid and the structure 
of the dilations can be exploited to obtain very sharp results. For irregular grids 
and unstructured dilations or if dilations are replaced by other transformations 
the situation is more complex and requires different techniques. One method is to 
use the regular case and try to obtain perturbations of the grid that preserve the 
frame structure. Another possibility is to obtain irregular samples of the continuous 
transform, that have the required properties. 

In this article we study frame decompositions of the space L 2 (R d ) using translations 
of a family of functions on irregular grids, and arbitrary dilations, and we even 
replace dilations by other transformations. 

Our approach is different and very general, allowing quite general constructions. We 
prove the existence of smooth time-frequency frame atoms in several variables. The 
setting includes as particular cases, wavelet frames on irregular lattices and with 
a set of dilations or transformations that do not have a group structure. Another 
particular case are non-harmonic Gabor frames with non-uniform covering of the 
Euclidean space. It also leads to new constructions of wavelet and Gabor frames 
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with regular lattice translates. One of the nice features of the proposed method is 
that it unifies different atomic decompositions. 

For the case of regular lattices Guido Weiss and his group |HDW02llHDW03llLab"02] 
developed a very fundamental program to characterize a large class of decomposi- 
tions of L 2 (R d ) through certain equations that the generators must satisfy. This is 
an important attempt to unify Gabor and wavelets decompositions. Other funda- 
mental construction of MRA wavelet frames on regular lattices can also be found 
in |CHS| . |CHSS03] . [CSOO . Our methods can be used to produce a substantial 
part of these systems. 

A set Q C R d is a wavelet set if the inverse Fourier transform of the characteristic 
function of the set is a wavelet. Wavelet sets, frame wavelet sets and methods 
for constructing such sets have been studied recently |BMM99] . |BL99] , [BLOlj . 
[BS03] . |DLS97| . |DLS98j . [HLOO] |Ola03j . [OS03] . Our methods give constructions 
of wavelet sets with translations on irregular grids. 

Let J and K be countable index sets. We consider families of functions {gj}jeJ C 
L 2 (R d ) and discrete sets X = {x j: k ■ j € J, k € K} C R d such that the collection 
{gj(x — xj,k) ■ j € J, k £ K} form a frame for L 2 (R d ). The wavelet case is obtained 
when gj = \ det(A)\^ 2 g o with A an expansive matrix and g a fixed atom. We 
want to stress here that our constructions are much more general, allowing for 
example a different invertible (not necessarily expansive) matrix Aj for each j G J. 
For the case of orthogonal wavelets, Yang Wang [Wan02 has recently considered 
wavelet sets associated with arbitrary families of invertible matrices and irregular 
sets of translates. He gave conditions for the existence of such wavelet sets and 
related them to spectral pairs. 

Irregular wavelet and Gabor frames also have been studied as perturbations of uni- 
form (lattice translate) frames and also as sampling of the continuous wavelet/ Ga- 
bor tra n sform. S ee |Bal97| , |BOHL03|, |Ohr96|, prr97]. |OFZf)l|, |CH97j, [GDH99] , 
[FZ95] , |FG89] , |Gro91| . |Gro93| . |HET)3] . [0592], [R595] , [SZOT)] . [SZOT] . [SZT)2j . 
[SZ03], [SZ03]. 

The approach in this article can be considered in the spirit of the classic construction 
in 1 dimension of smooth regular tight frames done by Daubechies, Grossmann 
and Meyer in [DGM86 . They found, for the case of uniform lattices, general 
conditions on a compactly supported smooth function h, in order that it generates 
a tight Gabor frame of L 2 (R) . In the affine case they found necessary and sufficient 
conditions for a band limited function in order that it generates a smooth wavelet 
frame. See also |HW89j . 

There were other related attempts to obtain atomic decompositions of functional 
spaces using very general systems. See for example FG85], |Fei87j in the context 
of locally compact groups. 

This paper is organized as follows: Section [2] introduces the notation and some 
preliminaries. Section [3] presents a Theorem on wavelet construction on arbitrary, 
sufficiently dense, but otherwise irregular grids and with arbitrary dilation or even 
invertible transformation matrices. Specific constructions of such wavelets are ob- 
tained in Section [4] first in the 1-D case and then in the multidimensional case. A 
general theory of frame atomic decomposition of L 2 (R d ) is obtained in Section pi 
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Using the concept of outer frame, reconstruction formulas for these atomic decom- 
positions are obtained in Section [6] 

2. Notation 

Throughout the paper J and K will denote countable index sets, and e x will stand 
for the function e^^) = e~~ l27rx '^. We will use p(E) to denote the Lebesgue measure 
of a measurable set E. 

A set Ti :— {hj}j£j of measurable functions on R d is called a Riesz partition of 
unity (RPU), if there exist constants < p < P < +oo such that 

(2.1) P <^2\hj(x)\ 2 < P a.e. x E M. d . 

Let iS = {Sj}j£j be a family of measurable subsets of K d . A Riesz partition of 
unity associated to S, is a set H := {hj}j^j of measurable functions, such that 

(1) Supp hj C Sj 

(2) There exist constants < p < P < +oo such that 

(2.2) P <^2\h 3 (x)\ 2 < P a.e. xeUjSj. 

jeJ 

Remarks. 

• If H = {hj} is a RPU, then H = {hj} is also a RPU. 

• If p = P = 1, we will say that H = {hj} is a regular partition of unity. 

• If the sets in S are essentially disjoint (i.e. n(Si (1 Sj) = 0,Vi ^ j), the 
family {hj — xs, } will yield a regular partition of unity associated to iS. 

• Every RPU TL = {hj} can be normalized to obtain a regular partition of 
unity by considering 

hj = (E,-N 2 ) 1/2 ' 

• Given a family S — {Sj}j^j of measurable sets on K d , define 

(2.3) p s (x) - #({j e J : x e Sj}) = ^X5 3 (^), 

jeJ 

where is the cardinal of the set B. The value ps = \\ps\\oo is called 

the covering index of S. 

We now recall the definition of frame for a given close subspace F of L 2 (R d ) . 

Definition 2.1. A set of functions {gj}je,J 1S a frame for F if gj £ F and there 
exist constants < m, M < +oo, such that 

(2.4) m||/|| 2 <^|</, ft >| 2 <Af||/|| 2 , V fGF. 

jeJ 
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For a measurable set Q C R we will denote by ICq the functions that have support 
in Q, and by Bq the functions whose Fourier transform has support in Q , i.e. 

(2.5) K Q := {/ e L 2 (M d ) : Supp(/) C Q} 

(2.6) £ Q := {/ e L 2 (R d ) : Supp(/) C Q}. 

We will need also the following definition. 

Definition 2.2. Let H be a Hilbert space. A collection of vectors {gj € H}j e j 
is an outer frame for a closed subspace -F of ff, if {-PF(5j)}jeJ is a frame for F 1 , 
where Pp is the orthogonal projection onto F, or equivalently, there exist constants 
< m, M < +oo, such that 

(2.7) m||/|| 2 <^|</, ff3 >| 2 <Af||/|| 2 , V fEF. 

Related definitions to the concept of outer frames appear in [L098 , [FW01 . 

Remark. Throughout the paper, we will use the following immediate and very useful 
fact about frames: 

If {gj}je.J is a frame for ICq, and V C Q, then {gj}jeJ is an outer frame for /Cy. 



3. Wavelets on arbitrary irregular grids and with arbitrary 
dilation matrices and other transformations 

Our first results concerns the construction of wavelets {| detAj\ 1 ^ 2 ip(AjX — Xj^) '■ 
j G J, k € K} with translates on the arbitrary irregular grid X = {xj.k G M d : 
k G K, j G J} and with an arbitrary countable family of invertible d x d matrices 
{Aj G GL d (R) : j G J}. 

Theorem 3.1 (Wavelets). Let Q C R d be a set of finite measure, h a function in 
L 2 (R d ) and A = {Aj G GLd(R) : j G J} a family of invertible matrices. 

For each j G J se£ Bj = (Aj) -1 , Sj = B^ 1 ^ = AjQ, h 3 = h{B r ) and let 

s = {s j ,jeJ}. 

Assume that S is a covering of M. d , TL is a RPU with bounds p and P and that 
Supp(/i) C Q. 

Consider X — {xa^ G R d : j G J, k G K} such that for each j G J, the set 
{e Xj k XQ '■ k G K} forms a frame for ICq with lower and upper frame bounds rrij 
and Mj respectively. If m := infj rrij > and M := sup^ Mj < +oo, then the 
collection 

{| det A j \ 1/2 ip(A J x - x jtk ) : j G J, fc G K} 

is a wavelet frame of L 2 (M. d ) with bounds mp and MP, generated by a single func- 
tion ip, where ip is the inverse Fourier transform of h. 

Proof. Since for each j G J we have that {e Xj k XQ '■ k G K} forms a frame for 
ICq with lower and upper frame bounds rrij and Mj respectively, an application of 
Part 2 of Corollary 5.2 for the matrix BJ 1 * shows that {\Bj\ 1 ^ 2 e Xj k (Bjw)xQ(Bju) 



k G K} forms a frame of K.$ with the same bounds. From the definition of Sj, 
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{([J>(Sj)) < " 1 ^ 2 ^ e (AT 1 Xj fc )( w )X5j ( w ) : k € K} is then a frame for JCsj with frame 
bounds mj[i(Q)~ 1 and Mj/i(Q) _1 . 

On the oth er si de, if Supp(/i) C Q then 7i is associated to S. So, we can apply 
Proposition 5.7 to 5,W and {^(S,-) - ' 1 ' 2 ) e^-i x . fc )XSj : & € J € ^; to conclude 



that, {fi(Sj)^ 2 hje^ A -i x , ^ : k G if,j G J} forms a frame of i 2 (IR d ) with lower 
frame bound mp/i(0;) _1 and upper frame bound M P^l(Q)^ 1 . This gives that, 

{\Bj\^ 2 h(BjUj)e Xjik {BjUj) : k G ifj G J} 

forms a frame of L 2 (M d ) with frame bounds mp (or mc) and MP. The theorem 
now follows from an application of the inverse Fourier transform. 

□ 

Remarks. 

(1) The set of matrices {Aj G GL^M.) : j G J} can be arbitrary and need not 
have a group structure. 

(2) The set {Aj G GL^M) : j G J} can also be chosen to have a simple 
structure. For example, J = Z 2 , ^4(i,j) = ii^-D- 7 where i? is a rotation and 
D a dilation matrix, will be used to construct directional wavelets. An even 
simpler example is J = Z, Aj — A J , where A is an invertible matrix which 
gives a construction of wavelet frames on M. d . 

(3) Note that h does not need to be compactly supported. 

We will use the theorem above to construct specific examples of wavelets, e.g., 
directional wavelets, isotropic wavelets, etc. 



Interesting particular cases of Theorem |3.1[ 

(1) x.j. k = x k y j g j. 

Let X = {xk G M. d : k G K} be such that {e Xk \Q,k G K} is a frame for 
K.q with frame bounds m and M. Then 

{IdetAj^^^AjX-Xk) : j G J,keK} 

forms a wavelet frame of L 2 (M. d ) with bounds mp and MP. 

(2) Aj = A^V j G J, with A G Gi d (M). 

Each of the following sets are wavelet frames of L 2 (JBL d ) with bounds mp 
and MP: 

{| det A\ j/2 tp(A j x - Xj, fe ) : j G J, fc G if }, and 
{| det A| J ' /2 V>(A J ':r - x k ) : j G J, k G if}. 

(1) If the set S is a tiling of M. d then the wavelets constructed above are 
Shannon-like wavelets, thus not well localized in space. To obtain well 
localized space-frequency wavelets, Tl must be constructed to be a smooth 
partition of unity, e.g., at least C 1 (K' i ) as demonstrated in the examples in 
Section |4j below. 
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(2) Reconstruction formulas for such wavelet frames are developed in Section 



4. Examples of wavelet frames on irregular lattices and with 
arbitrary set of dilation matrices and other transformations 



4.1. Density and Separation. To be able to use Theorem |3.1| to construct con- 
crete examples of wavelet frames on irregular grids, we first need to construct ex- 
ponential frames {&x k XD} (also called Fourier frames) for ICd- Exponential frames 
play a central role in sampling theory for Paley- Wiener spaces (also known as spaces 
of band-limited functions) . 

The density of a set X = {xk C M. d : k G K} and separateness of the points in X 
play a fundamental role for finding exponential frames {e Xk xo} for ICd- 

Definition 4.1. A sequence X = {xk ■ k € K} is separated if 

inf \\x k - xi\\ > 0. 

k^l 

There are many notions for the density of a set X. We start with three definitions 
that are due to Beurling. 

Definition 4.2. 

(1) A lower uniform density D~(X) of a separated sequence X C R d is defined 
as 

v-(r) 

r — >oo 



D~{X) = lim 



(2r) d 

where v~(r) := min # (X n (y + [-r, r] d ) , where #(Z) denotes the cardi- 
nal of the set Z. 

(2) An upper uniform density D + (X) of a separated sequence X is defined as 

D+(X) = lim 

r— voo (2r) 

where z/ + (r) := max# (X fl (y + [— r, r] d ) . 

(3) If D-(X) = D+(X) = D(X), then X is said to have uniform Beurling 
density D(X). 

Remark. The limits in the definitions of D~(X) and D + (X) exist (see [BW99 ). 

As an example, let X C K be separated and assume that there exists L > such 
that \x k - || < L, for all k e Z. Then D-(X) = D+(X) = d. For the one 
dimensional case, Beurling proved the following Theorem. 

Theorem 4.3. (Beurling) Let X C K be separated, a > and = [— |, §]. // 
a < _D - (X) f/iera {e a;( .xn} is a frame for /Co- 



This previous result however is only valid in one dimension. For higher dimensions, 
Beurling introduced the following notion: 
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Definition 4.4. The gap p of the set X — {x k : k G K} is defined as 
9 



p[X) = inf J r > : |J B r (.T fe ) = R d I 

I k£K J 



Equivalently, the gap p can be defined as 

p = p(X) = sup inf \x — Xk\. 

x eR dXk€X 

It is not difficult to show that if X has gap p, then D~(X) > i. For a separated 
set X, and for the case where f2 is the ball B r (0) of radius r centered at the origin, 
Beurling [Beu66 proved the following result: 

Theorem 4.5 (Beurling). Let X C R d be separated, and = B r (0). If rp < 1/4, 
then {e Xk xn} is a frame for ICq. 

For a very clear exposition of some of the Beurling density results see |BW99| . 
4.2. Wavelet frames in 1-D. The construction in the following theorem (which 



is a particular case of Theorem 3.1 ), generalizes a similar result of [DGM86 to the 



irregular case. See also |Gro93j . 

Theorem 4.6. Let Q = [-1,-1/2] U [1/2,1], < e < 1/2, and let h + be a real 
valued function such that Q\ := Supp h+ C [i — e, 1 + e], \h+\ < 1, and < c < \h+\ 
on [1/2, 1]. Assume that for each j G Z, the sequence Xj — {%j,k} kG % is separated 
and that D~(Xj) > 2 J+1 (1 + e). Then for each j, the set {e 2 i Xj k XQ e '■ x j.k € 
Xj : k G Z} is a frame of K,q c , where Q e := Q\ U (— Q\). If furthermore the sets 
Xj = {xj,k} kG % are chosen such that the frame bounds rrij and Mj satisfy infj rrij — 
m > and supj Mj = M < +oo, then the set {2 j / 2 tp(2 j (■ - x ]M )) : j G Z, k G Z} 
where ip(x) — 2Re{h + {x)) is a wavelet frame for L 2 (M). 

Remark. The wavelet frame constructed in the theorem above is of the form 
{2 J '/ 2| 0(2 J '(- — Xj t k)) '■ j G Z, fc G Z} which is slightly different form than the one 
constructed in Theorem |3.1| This discrepancy is due to a convenient choice of the 
irregular set Xj* — {xj t k} that we have adopted in the statement of the theorem 
above. 

Note that the wavelets constructed in the theorem above are real and symmetric. 
Actually, if one wants tp with good decay, h + can be easily constructed to be C r , 
r > 1, even C°°. 

As a corollary, if we choose the sampling sets Xj to be nested, i.e., Xj C Xj + i, we 
get 



Corollary 4.7. Let Q, e, and h + be as in Theorem 4-6 Assume that the sequences 
Xj = {£j,fc}fc<=Z are separated and such that Xj t k — Xj+i.ik- If D~(Xq) > 2(1 + e), 
then for each j , the set {e2j Xj k XQ e '■ x j,k *= Xj} is a frame for K-Q e . If furthermore 
the frame bounds rrij and Mj satisfy infj rrij = to > and sup^ Mj — M < oo, then 
the set {2 3 '' 2 ip(2r' (■ — Xj^)) '■ j G Z, k G Z} where ip(x) — 2Re(h + (x)) is a wavelet 
frame for L 2 (E) . 
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Proof. Since Xj^ = Xj+i,2k we have that #(Xj+i H [—r, r]) > 2#(X,- n [—r, r] — 1). 
Thus D-(X j+i ) > 2D-(Xj). But D~(X ) > 2(l+e), therefore £>-(X,-) > 2 j+1 (l+ 



e). The corollary then follows directly from Theorem 4.6 □ 



From Theorems |4.6| and |3.1[ we immediately get the following Corollary. 



Corollary 4.8. Let Q, e, and h + be as in Theorem 4-6 Assume that the set X = 
{^fejfegZ * s se P ara t e d and that D~{X) > 2(1 + e). Then, the set {e Xk XQ € ■ Xk G X} 
is a frame of JCq c , and the set of functions \2?l 2 i\)(2? ■ —x^) : j G Z, k G Z} where 
ip(x) = 2Re(h+(x)) is a wavelet frame for L 2 (M). 



4.2.1. Examples. 



4.7 



with e = and h 4 



(1) Shannon-type wavelet frames : We use Corollary 
X[i/2,i] to get wavelet frames of the form 

i>x j>h ,j = {2 _j/2 cos(2 _: '~ :L 37r(a; - sinc(2 _: '~ 1 7r(a: - Xj,k))}- 

These wavelets are not well localized since the decay at oo is O(|o;| _1 ). 

(2) Shannon-type wavelet bases : If we choose Xj such that \xj k — 2 J k\ < Lj < 

V k G Z, then by Kadec's 1/4-Theorem, we immediately get that 
iP Xjkd = {2-J/ 2 cos(2-J- 1 37r(a;-a; J - fc ))sinc(2-J- 1 7r(x-a; J , fc )) : j G Z, fee 
Z} constructed above form a wavelet Riesz basis for L 2 (M). 

(3) Well localized wavelet frames: For faster decay of the wavelet frames, we 
choose h + to be a smoother function. Let f3 n = X[o,i] * • ■ • * X[o,i] be the 
B-spline of degree n (note that Supp/3„ = [0,n+ 1]). Let e = 1/4, and 
/i+(^) = n(3 n -\ ((£, — l/4)n)). Then we get a wavelet frame of the form 

4> Xj>kd = {2^/ 2 cos(2^'- 1 3^( a; - ^- fc ))sinc"(2-^(x - x j<k ))}. 
For this case the wavelets decay as 0(|x|~"). 

4.3. Examples of wavelet frames in K d . 

(1) For W 2 , let X = {xk '■ k G Z}, and let Y = {y t : I G Z}. If D~(Jf ) > 2 and 



D (Y) > 2, then using Proposition 5.6 below for product frames, the set 
{e( Xk ,yi) '■ (k,l) G Z 2 } form a frame for /C[_ 11 ]2. Let Q := {(x,y) G K 2 : 
1/2 < x 2 + y 2 < 1}, and A = 21, then R 2 = UjA j Q. We can then use 
to construct wavelet frames for L 2 (R 2 ): 



Theorem 



3.1 



• Shannon-type radial wavelets: Let h = \q, then h is radial. Thus 
the function ip defined as ip = h satisfies tj){x,y) — g(r), where r = 
(x 2 + y 2 ) 1 / 2 . We then cons truct the Shannon-like wavelet frame for 



L 2 (R 2 ), as in Theorem 3.1 A related construction of non-separable 
radial Shannon-type frame wavelets and multiwavelets can be found 
in jPCKKHlj . and |PCKKH2j . 

Well localized radial wavelets: To construct wavelet frames with poly- 
nomial decay in space, we let h^i,^) — n/3 n - i ((£ i + Cl — 1 /4)n) , 



and construct the wavelet frames using Theorem 3.1 (see Figure [T]) 
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Figure 1. Radial wavelet frames that are well localized in space. 



The points {(xk, yi) '■ k, I £ Z 2 } lie on an irregular grid of the form X x Y. 
However, we may be interested in points Z = {(xk,yi) CR : (fc,Z) £ Z } 
that do not lie on irregular grids of the form X x Y . For this case, the 
same constructions above can be used to form wavelets frame for L 2 (K 2 ), 
as long as the gap p(X) < 

As in the 1-D examples above, we can also use Corollary 4.7 to construct 
wavelets on irregular grids satisfying Xj C Xj + i. 

Directional wavelet frames: We can easily construct directional wavelet 
frames as follows: Let Q\ be a region defined by Qi — {(x, y) £ M 2 : x = 
rcos(6),y = rsin(0),l/2 < r < 1, |0| < f }, and define Q = (-Qi) U Q x . 
Let A — 21, and R be the matrix of a rotation by an angle ir/4. Let ip be 
such that ^ = xq, then we obtain the wavelet frame for L 2 (M 2 ) of the form 
= 4-JV 2 ^(2- J ' 1 ^ 2 • -x Jtk ) : j = (ji,j 2 ) eZx {0,1,2,3},* £ Z}. 
The index j\ codes for the resolution of the wavelet, while the index codes 
for four possible directions. Thus the wavelet frame coefficients encode time 
scale as well as directional information. Clearly one can choose any number 
of directions and adapt the previous construction. An obvious modification 
as shown in Figure [2j yields wavelet frames with polynomial decay. Very 
nice constructions of smooth directional wavelet frames on regular grids 
were obtained before in |AHNV01| lADH+03] 



Spiral In this example we will define a dilation covering by spiral annulus 
sectors. 

Let a,b > 1, and T the spiral curve defined by 

r(t) = (a* cos(bt), a* sin(bt)) t £ E. 



For a £ K define R a to be the rotation of angle a : R a = [ c °^2 



- sin a 
cos a . 

The curve T satisfies: 



T(t + a) = a a R ba T(t). 

Note that for positive a the matrix A = a a Rb a is expansive. 

Now we are ready to define the covering elements. Set b — 2ir and a = 
for some integer to > 2 so that A m = aid- Define the spiral annulus sector 
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Figure 2. Well- localized directional wavelet: The regions Q = 
Q 1 U (— Q 1 ) and Q t — Q\ U {—Q\) that can be used to construct 
well localized directional wavelets. 




Figure 3. Spiral wavelet frames. 

Q = {x e K 2 : x = XT(/3), l<A<a,0</3<a} (see Figure ^. So Q is 
compact and {A>Q : j G Z} is a disjoint covering of M d \{0}. 

Choose e > sufficiently small and h a smooth function that does not 
vanish in Q and with support in Q e . Define ip = h. Select a separated set 

X = i X k} keZ C M 2 SUCh that < 2 diam(Q E ) 1 

The set {a^ m ip(a j / m R-2„j/ m (x - x k ),k G Z,j G Z} form a wavelet 
frame of L 2 (K 2 ) generated by a single wavelet -0 that is band-limited, with 
good decay and directional in frequency. 

(6) Obviously, all the constructions above can be generalized to Mr for any 
dimension d > 2. 

Remark. Some of the wavelet frames may be associated with MRAs. For exam- 
ple, the so called Shannon wavelet frame constructed above is associated with the 
Shannon MRA V 3 = {/ G L 2 (R) : Supp/ C [-2" 5 '- 1 , 2~ j - 1 ]}, j G Z. In general 
however, the precise relation needs further investigation. 
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5. General results on the construction of time-frequency atoms 

In this section we will develop a general method for constructing time-frequency 
frame atoms in several variables. This construction allows us to construct the 
previously introduced wavelet frames on irregular grids and with arbitrary dilation 
matrices or other types of transformations. It also allows us to construct non- 
harmonic Gabor frames on non-uniform coverings of R d as described in Section 5.] 
below. 

Let U and V be non-empty open subsets of M. d , T : U — ► V is an invertible C 1 
map, with C 1 inverse S : V — * U, i.e., T is a C 1 homcomorphism with C 1 inverse 
S : V -> U. Define a := \nf vBU \detT'(y)\- 1 and f3 := sup„ etr | detT'iy)]' 1 , where 



T(y) denotes the derivative of T at y, T'(y) = 



Proposition 



We have the following 



Proposition 5.1. With the above notation, assume that a is positive and (3 is finite, 
then if {gj}je.J * s a frame for K,y with frame bounds m and M , then {gj o T} je j 
is a frame for K,jj with frame bounds am and (3M . 

Proof. Assume that / 6 K,\j- Then for j £ J, by an application of the change of 
variables formula 

(5.1) <f, gj oT>= [ fgJ^T = [ (foS)\detS'\gJ =<(foS)\detS'\, gj > . 

Ju Jv 

Since | det S'\ is finite and bounded away from zero, (/ o S) \ det S'\ is in JCv- 
Using that {gj}jej is a frame of /Cy with bounds m, M we have 

(5.2) ]T| < (foS)\detS'\, 9j > | 2 < M\\(f oS) \detS'\ \\ 2 . 

Now, applying again the change of variables theorem, we obtain 

||(/o5)|det5'||| 2 = / |/oS| 2 |detS'| 2 = / |/| 2 | det (5' o T)| 2 | detT'|. 
Jv Ju 

Since S = T _1 we have that det[(5' o T)(y)] = (det[T'(2/)]) _1 for all y£U. Thus 

(5.3) ||(/o5)|det5'||| 2 = [ l/Hdetrr^Afll/H 2 . 

Ju 



That is, from ( 5.1 ), ( pT2| and ( |5T3[ ) we have 



^| <f, 9j oT> \'<(3M\\f\\ 



This proves the upper inequality for {gj o T}j e j. The lower inequality is obtained 
in the same way, with the obvious modifications. □ 



Proposition 5.1 establishes that T defines an isomorphism IIt : K-u — > JCy defined 
by 11^(3) = g o T, thereby transforming frames into frames. By taking T to be a 
translation or a dilation, we get the well known result: 

Corollary 5.2 (Translation and dilation of frames). Let Q CK d be an open subset 
ofR d . Let y£R d be any point, and A £ GL^R) an invertible matrix. We have, 
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(1) Translation {gj}jej is a frame for ICq with bounds m and M, if and only 
if {gj(- — y)}j£j is a frame for 1Cq +v with the same bounds. 

(2) Dilation {gj}je.J is a frame for JCq with bounds m and M, if and only if 
{| det A\~ 1 / 2 gj(A~ 1 -)}j l zj is a frame for K-aq with the same bounds. 

Remark. In fact the corollary remains true if we only assume that Q is measurable. 

Proof. Part[l] is a direct application of the Proposition, for the case that T(x) — 
x - y. 

For Part [2J and the transformation T(x) — A~ x x, the Proposition tells us that 
{gj(A~ 1 -)}j e j is a frame with frame bounds |detA|m and |det^4|M, therefore 
dividing each function by | det A^ 1 / 2 we obtain the result. □ 

For the next theorem, we need to introduce some definitions. Let TL = {hj : j G J} 
be a RPU with bounds p and P. For each j £ J set 

(5.4) Wj = Closure^ihjf : f S L 2 (R d )}. 
Let < c < p, and define 

(5.5) Q 3 ; = Qj(c) = {x e R d : Ih^x)] 2 > c}, for each j £ J. 

For a given c, we discard all those j such that Qj has measure zero. Note that if 
Jo = {j £ J : m(Qj) > 0}' then we can only claim that TIq = {hjXQ^j G Jo} is a 
RPU associated to {Qj}je.J w ith constants c and P. 

Theorem 5.3. Let < c < p. 

(1) Assume that {gj_k}k<EK is a frame for Wj with lower and upper frame bounds 
riij and Mj respectively. If m :— infj mj > and M :— sup^ Mj < +oo. 
then {hjgj^k '■ j G J, k G A'} is a frame for L 2 (JH d ), with frame bounds pm 
and PM. 

(2) Assume that for each j G Jo, {gj,k}keK is a frame for ICqj with lower and 
upper frame bounds rrij and Mj respectively. If m := inf j € j rrij > and 
M := sup Jg j o Mj < +oo, then {hjgj,k : j G Jo , & G A"} is a frame for 
fcuQj, with frame bounds cm and PM. 

Proof. 

(1) Given / G L 2 (R d ) and denoting by f 3 =hjfe Wj we will first show that 

(5.6) p||/|| 2 < El^ll 2 ^ P ll/H 2 < 

j 

for 

pII/II 2 = J P \f\ 2 < /Ei^-i 2 i/i 2 = E / im 2 = Eh/.ii 2 ' 

3 j j 

where we used dominated convergence for the interchange of the integral 
with the sum. The other inequality is analogous. 
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For each j € J and / € J 2 (IR d ), we use the fact that {<7j,fc}fc is a frame 
for Wj, and that < f,hjgj^ >=< hjfi9j.k >, to obtain 



HIMI 2 < "b-ll^/ll 2 < E I < /.Ma* > I 



2 



< MjWhjff < M\\hjff. 

So summing over j 

pm\\f\\ 2 < mj2\m 2 < EEi</'^, fe >i 2 

j 3 k 

< Mj2\\fA 2 < PM\\f\\ 2 . 

j 

(2) For the second case we observe that the set A = {\j = hjXQ s ,j £ Jo} forms 
a RPU associated to {Qj,j € Jo} with bounds c and P, for if a; € UQj, 

c<c#({j6J :a;eQ 3 -})< E l^( x )| 2 

= Ei /i i( a: )i 2 ^^)<Ei^( a; )i 2 ^ p - 

ieJo ieJo 
Furthermore, by hypothesis, 

£q 3 - = {^/ : / e i 2 (M d )}, j e Jo. 



These subspaces JCq^ correspond to the subspaces Wj defined in (5.4) 
for the RPU A, and therefore we can use the previous result (applied to 
L 2 (UQj) instead of L 2 (R d )) to conclude that {Xjg^k ■ j G J,k G K} 
is a frame for JC\ J Q j — L 2 (UQj). The proof is complete by noting that 
Ajffj,*; = hjgj : k,j € J Q . 

□ 

Remarks. 

• Note that in the previous theorem, instead of choosing a frame for the 
subspaces Wj, we could have chosen any collection of functions of L 2 (M. d ) 
that form an outer frame for Wj . 

• If h is a bounded function and Q = Supp h, then it is easy to see that 
closure L 2{/i/ : / e L 2 (R d )} = JC Q if and only if ^(Q) = fi({x : \h(x)\ > 0}). 
So the spaces Wj defined in (5.4 1 will coincide in most of the cases with 

^Supp(/i J )- 

As a very important particular case of the previous theorem, we have the following 
Corollaries. 

Corollary 5.4. Let S = {Sj C R d : j S J} be a family of subsets of R d , not 
necessarily disjoint, and let TL — {hj} be a RPU with constants p, and P not 
necessarily associated to S . Assume that {gj,k}keK is a frame for ICsj with lower 
and upper frame bounds irij and Mj respectively. If m := inf j rrij > and M :— 
swpj Mj < +oo, then 
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(1) IfTtis associated to S (i.e. Supp hj C Sj,j G J), then {hjgj^k ■ j G J,k G 
if} is a frame for /C Us , frame bounds pm and PM. 

(2) If instead \hj(x)\ 2 > c Vie G S^,j G J, iften also {hjgj,k '■ j G J, k G if} «s 
a frame for JC u Sj, with frame bounds cm and PM. 

Proof. Denote by Pyy the orthogonal projection on the subspace Wj defined in 



(5.4 1. Since {gj,k '■ k G if} forms a frame for X^. with bounds mj and Mj, then 



{Pwi9j,k '■ k G if} forms a frame for Wj with the same bounds. So, by Theorem 5.3 
{hjPw Sj.k '■ j ' G J, k G if} forms a frame of /Cys with frame bounds mp and MP 
We obtain part (1) of the Corollary observing that: 

< /, (.'/.,./. ! >=< hjf,P Wj (g j>k ) >=< I'jf-Uj.i. >=< /• /',</.; /. > • 

The second part is a consequence of the fact that {hjXs, ■ j G J} is a RPU 
associated to 5 with bounds c and P and that Qj (c) = Sj . Therefore we can apply 
the second part of Theorem |5.3[ □ 

Remarks. 

• For the first case, the set {hjgj^k '■ k G if} is not necessarily a frame for 

or even for Wj C fCs, , even though {hjgj t k '■ j G J, k G if} is a frame 
for /C Us . } . 

• In the second case, the subspace Wj is contained in K,Sj- 

• Note that if Ti. is a regular partition of unity, then the frame bounds for 
the frame constructed in the Corollary, are m and M. 

Corollary 5.5. Let Ti = {hj : j G J} be a RPU with constants p, and P and 
{Sj C R d : j G J} and {Qj C R d : j G J}&e coverings of R d such that Qj C Sj 
and |/ij(a;)| 2 > c a.e. x G for all j € J and some constant c > 0. Assume 
that {gj.k}keK is a frame for K-Sj with lower and upper frame bounds mj and Mj 
respectively. If rn := inf j nij > and M :— sup j Mj < +oo, then {hjgj^k '■ j G 
J, k G if} is a frame for L 2 (JH d ), with frame bounds cm and PM . 

Proof. As in the proof of the Theorem, we have the inequalities 

(5.7) pII/II 2 < £llMI 2 ^ p ll/H 2 ' 

j 

On the other side, 

< /• ..//;.; >=< hjfi93,kXS s >=< l'.f\ - •.'/;./• >, 

and then 

EE I < /> h i9j,k > I 2 < Mj^Wh.fxsM < Mj2\\hjf\\l < MP\\f\\l 

j k j j 

and 

EE I < f> h i9i,k > I 2 > mE WhjfXBj Hi > mcY, \\fx Q) \\l > mc\\f\\l 

j k j j 

a 
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The following proposition is a direct application of Fubini's theorem, and allows 
us to construct frames in a product space. By induction the Theorem can be 
generalized to hold for any finite number of factors. 

Proposition 5.6 (Product of Frames). Let E\ and E2 be measurable subsets of 
K 1 and K 2 respectively, and let {hj}j e j and {gk}keK be frames for JCe 1 andK,E 2 
with frame bounds mi, Mi and m2, M 2 . Then {hjgk}j t keJxK is a frame of ICe 1 xe 2 
with frame bounds m — m\m,2, M — M\M2- 



Proof. For any function / in JCex x e 2 > 



EE 

j£,J keK 



f(x 1 ,x 2 )h j (x 1 )g k (x 2 )dx 1 dX' e 



E x xE 2 



EE 

jeJ \k€K 



E/2 \*J E± 



f(xi,x 2 )h j (xi)dx 1 gk{x 2 )dx2 



i E 2 j^j i Ex 



f(xi,x 2 )hj(xi)dxi 



dxo 



> m\m 2 



E 2 JE 



\f(xi,x 2 )\ 2 dxidx2, 



which yields the lower frame bound m = ra\m2- The upper frame bound M — 
M1M2 can be obtained in a similar fashion. □ 



5.1. Construction of time-frequency atoms on arbitrary irregular grids 
and with arbitrary dilation matrices and other transformations. We now 

particularize our previous results to frames of the form {he Xj h , j € J,k € K}, where 
h is a fixed function. Using the Fourier transform, these types of frames allow us 
to construct wavelets {\ det Aj\ 1 ' 2 ip(AjX — Xj y k) ■ j € J,k € K} with translates on 
the arbitrary irregular grid X and with an arbitrary countable family of invertible 



d x d matrices {Aj 6 GL £ ;(IR) : j £ J} (cf. Theorem 3.1|. First as a particular 



case of Corollary |5.4[ we obtain the following Proposition. 

Proposition 5.7. Assume that S = {Sj : j € J} forms a covering ofM. d , and let 
TL = {hj : j G J} be a RPU with bounds p and P associated to S. Assume also 
that {/x(S , J -) _1 ' 2 e x . k XSj '■ k € K} * s a frame for Kg. with lower and upper frame 
bounds rrij and Mj respectively. If m :— infj mj > and M :— sup^ Mj < +00, 
then 

{/'(* ; ) 1 : je J,keK} 

is a frame for L 2 (R d ) with frame bounds mp and MP. 

Remark. If xj t k, Sj, and hj are chosen such that Xj.k = ak, k G 7L d , a G M, 
Sj = S + j, j e T and hj = h(- + j), j G T, where T is a lattice in R d , then 
we obtain the standard Gabor or Weyl-Heisenberg frames. Thus in general, the 
construction above can be viewed as non-harmonic Gabor frames with variable 
windows hj. 

The following wavelet frame construction is a direct application of the previous 
proposition. Taking ipj = hj, the set {ipj tXj k — — Xj^) ■ j G J, k G K} is a 
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wavelet frame for L 2 (M. d ) with frame bounds mp and MP. Therefore, there exists 
a dual frame tpj, Xj k such that 

/ = EE < /> > <<w v ^ e i2 ( M ")- 



We can also use Proposition |5.7| and Beurling Theorem |4.3| to obtain a non- 
harmonic Gabor frame as the following example shows: 

Example. [Non-harmonic Gabor frames] Let Q — [1, 3] and Let — X[o,i] * 
X[o,i] * X[o,i] * X[Q,i] be the B-splinc of degree 3. Clearly, K = l) je %Q + j. Let 
Xj = {£j,fc} feg 2 be se ^ s m ^ chosen such that for each j, D~(Xj) > 2 and such that 
the frame bounds rrij and Mj satisfy infj mj = m > and sup^ Mj — M < +oo. 



Then using Proposition |5.7| and Beurling Theorem 4.3 we obtain a non-harmonic 
Gabor frame of the form {P3(t—j)e~ l27rXj ' kUJ : j S Z, fc € Z}. Obviously, we can also 
use a non-uniform partition of K and get generalized non-harmonic Gabor frames 
as discussed earlier. 

Remark. If Q C M d be a measurable subset of W 1 , and if the family of sets S is 
defined by means of expanding or contracting Q, then we obtain the theorems of 
Section [3] as particular cases of Proposition |5.7| 



5.2. Construction of exponential frames. By Corollaries |5.2| and |5.4| we can 

construct a frame for Kq C R d starting fr om a frame for K. d where D is any subset 
of M. d with nonempty interior. Hence to build a frame for ICq it is enough to start 
with a frame for K-u, where U is an open disk. Specifically, since any bounded 
measurable set can be covered by a finite number of translates of U, we can use 



Corollary 5.2 1) to find a frame for ICq. We can also expand U until it covers Q 



and use Corollary 5.2 '2). This shows that there are many ways to construct frames 
for ICq starting from a frame for ICu- Obviously, the particular construction will 
depend on the application. 

We will describe now two particular constructions. 

(1) Given a frame for Kd where D is a measurable subset of K d , we will 
construct a frame for Kq. Let T be a regular lattice in R d (i.e. L = RZ d , 
where R is an invertible d x d matrix with real entries), and let D be a 
measurable subset of R d such that K d = |J (D + 7) with a finite covering 

index. Let Q be a measurable subset of M. d and define Q 1 :— Q (1 (D + 7) 



for 7 e T. Let A := {7 e T : /i(Q 7 ) > 0}. By Corollary pT2| if {e Xk XD ■ 
k € K} is a frame for TCo, then {e Xk XD+-y '■ k S K} is a frame for JCfj +7 . 
Hence, {e Xk XQ^ ■ k £ K} is also a frame for ICq . Therefore, by Corollary 
5.4 {e Xk XQ~, '■ 7 £ A,fc S if} is a frame for ICq. As an example, when 
Q is a measurable subset of R , I? = [0, l] d , and L = Z d , we have that if 
{e Xk X[a,i] d : fc e if } is a frame for L 2 ([0,l] d ), then {e Xfe XQ 7 : 7 e A, fc G 
if} is a frame for Kq (recall that Q 7 = Q n ([0, l] d + 7)). 

It is easy to see that the covering requirement R d = [J (D + 7) in 

the previous construction is not restrictive. Specifically, we only need 
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7), where a is any positive real number. Furthermore, 



R d = U (aD 

■yer 

the construction remains valid if for each Q 7 we choose a different set 
{x 7 k} such that E 1 := {e x k \D '■ 7 G A, fe S if} is a frame for /C^, and 
m := inf 7 m 7 > and M :— sup 7 M 7 < +00, where m 7 and M 7 are the 
lower and upper frame bounds of E 1 . Finally, note that if Q = M. d we 
obtain a frame for L 2 (R d ). 
(2) If Q is bounded we can construct frames for Kq using Theorem 4.5 Let 
S = diam(Q), and x Q e R d such that Q C B(a3 , (5). Let X = Izfc, k e K } 
be such that p(X) < i. Then using Beurlings Theorem (4.5 1 we obtain 
that {ea; fc XB(o,<5)i ^ £ K} is a frame of K B ( s)- So, using Corollary 5.2 1), 
{ e x k XB(x ,&)->k e K} is an outer frame of Kq, and {e Xk XQ,k e K} is a 
frame for i^Q. 



5.3. Existence and construction of Riesz partitions of unity. In view of 
the previous results, we will be interested in constructing particular kinds of Riesz 
partitions of unity associated to special coverings of the space. The next results 
provide the necessary tools to accomplish this task. 

If A is a d x d matrix, we will say that A is expansive, if |A| > 1 for every eigenvalue 
A of A. 

We will use the following known result (see for example [HJ91J,pg. 297). 

Lemma 5.8. Let B be in <C dxd and e > 0. There exists a matrix norm |||'||| such 
that 

s(B) < \\\B\\\ < s(B) + e, 
where s(B) is the spectral radius of the matrix B, and there exists a norm \\ ■ | in 
C d such that 

\\Bx\\<\\\B\\\ \\x\\. 

As a consequence of this lemma, if A is an expansive matrix, then there exists a 
norm || • || in C d such that 

IIA-^H < c\\x\\ 0<c<l, 

and therefore 

\\Ax\\ > c'\\x\\ c'>l. 

In particular for every x G C d , 

(5.8) lim = and if x ^ lim \\A 3 x\\ = +00. 

Proposition 5.9. Let A be a d x d expansive matrix and V C R d a bounded set 
such that 

(i) there exists e > such that B(0,e) n V — 0. 

(ii) U ]eZ AW = R d \{0}. 

Then pa,v, the covering index of the family {A J V}j eZ , is finite, i.e. there exists 
an integer n > 1 such that 1 < pa,v ^ n - 
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Proof. For x £ M. d define 



z, 



+ _ 



{j > : x £ AV) and z~ = {j < : x £ A j V}. 



Since pa,v < ess sup,,. (#(2+)) + ess sup a ,(#(z~)), it is enough to prove that #(z+) 
and #(z~) are uniformly bounded in R d . We will see that 4^{ z t) is uniformly 
bounded. A similar argument proves the claim for 

It is easy to see that 

#(z+) = #(z+ J VseZ, V,eM d . 

Thus, if #(z+) is bounded on V, then, by ([h| #0+) is bounded in M d \ {0} with 
the same bound. Now, using we see that there exist < c\ < c 2 < +oo such 
that ci < < C2, for all x £ V. Fix a; € V. If j £ z+, then there exists i> € V 
such that x = A^v and 

C2 > n = \\a^\\ > > " 



\\A~q ~ \\A-if 

and therefore ||^1~ ; '| > f^- But since ||^4"- ; || — > 0, necessarily there exists jo such 
that for every s > jo, ||^4~ s || < Hence for every x £ V, #(z£) < jo- O 

The next proposition shows the construction of a large class of Riesz partitions of 
unity, for families of sets obtained by dilation of a compact set. 

Proposition 5.10. Let Q C R be a compact set and A a d x d expansive matrix 
such that 

(i) Q 

(a) U je z Aj Q = ^ d \{°}- 

Let h be any measurable function, and < C\ < c 2 < +oo some constants such that 

(a) < \h\ 2 < c 2 

(b) < ci < |fr| 2 on Q 

fa) /i = on M d \ Q e , where < e < d(0, Q) and Q e = {x £ R d : d(x, Q) < e}. 

Then the family of functions {hj(-) — h(A~ J * s a associated to {A J Q e }. 

Proof. If x ^ 0, by ([h]) there exists j £ Z and q £ Q such that = a;, so by |b|, 

ci < |«(g)| 2 = |/i(A- J 'a;)| 2 . Thus |n(A- s x)| 2 > c x . 

seZ 

Now by Proposition |K9j the covering index pa,q c of the family {A 3 Q e } is finite. 
Using Q and (JTTJ> , we see that Supp(/ij) C A^Q e and since < |ft.j(x)| 2 < c 2 V x, 
we obtain that 

\ h j( x )\ 2 ^ PA,Q e C 2 , 

jeZ 

which proves the proposition. □ 

Remark. This proposition generalizes easily to the case where we replace for each 
j, A } by an invertible matrix Aj in such a way that {AjQ} .^g is a covering of 

R d \ {0} with finite index. 
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The next Lemma shows, that the assumption of having a compact set that covers 
M. d by dilations is actually necessary, if one wants to construct a RPU associated 
to a covering of M. d of the form S = {A^V, j G Z}, where V is a bounded open set. 

Lemma 5.11. If A is invertible, and V is a bounded open set such that 

(1 ) there exists e > with B(0, e) H V = 

(2) {J jeZ AW = M d \{0} 

then there exists a compact set Q C V such that {J .^A^Q D R d \ {0}. 

Proof. It is clearly enough to find Q such that UjA^Q D V. Let 

Q n = \ x e V : d(x,dV) > - 
I n 

We will prove (by the contradiction) that for some n > 1 , Q n covers V by dilations 

by A. 

Assume that for each n > 1, there exists x n £ V such that x n £ UjA^Qn. Since V 
is compact, and d(x n ,dV) < — , there exists a subsequence {x nk } and x € dV such 
that x nk — > x. By and ([2| there exists jo such that x € A^°V, and since A^ ^ 
is open, x„ fc € A 70 ^ for k > k . Set y = A~i°x and yjt = A^°x nk , then y,y/c € V 
for fc > fc . Choose e small enough such that B(y,e) C V, and so d(y,dV) > e. 
Let mo be such that if to > too, then y m G -B(y, |). 

Then for z E B(y, |), and v £ dV we have 

e - d(y, z) < cZ(y, u) - d(y, z) < d(z, v), 

and thus 

£ -<d{z,dV). 

So -B(y, |) C Q n for all n such that — < § . This contradicts our assumption that 
x « m = ^ Jo l/n m ^ Qn mI since, if to > to and ~ < § , we have y m S i?(y, |) C Q„ m , 
and so x n € A? Q„ m . □ 



We show next that it is not difficult to construct bounded sets that cover K d by 
dilations. 

Lemma 5.12. Let V C R d be a bounded set such that £ V° , and A an expansive 
d x d matrix. Set Q — AV \ V , then {A^Q,j £ Z} is a covering o/K d \ {0} with 
finite covering index. Furthermore, ifVC AV then the sets {A J Q} are disjoint. 

Proof. Choose e > such that £?(0, s) C V. Let x £ R d \ {0} be an arbitrary point. 
By equation (5.8 1 lim^oo ||A _:, x|| = 0, for some norm || • || in C d . 

So there exists a positive integer n such that Vj > n, y :— A~^x £ 5(0, e), i.e., 
x £ A j V,Vj > n. 

Since V is bounded, and ||x|| > 0, there exists jo S ^ such that x £ A : >" +1 V and 
.x £ A 3 V, Vj < jo, i-e., x € A jo+1 V \ A jo V = A^Q. 



The finiteness of the covering index follows from Proposition |5.9| and if V C AV 
the disjointness property follows immediately. □ 
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Remark. 

• If we want to obtain coverings of R d \ {0} with compact sets, we can choose 
the set Q where Q is the set constructed in the previous Lemma. Clearly the 
sets {A>Q : j e J} cover R d \ {0}, and if fx(dV) = 0, then {A>Q : j E J} 
are almost disjoint, i.e. ^(A^Q n A k Q) = 0, j ^ k. 

• Even though for any set V it is true that A j (AV \ V) n A> +1 (AV \ V) = 0, 
in general it is not true that the family {j4 j '(AV \ V)} is disjoint, as the 
following example shows: Let A = \% ~q 2 ] and V the rectangle [—3,3] x 
[-1,1]. 



5.4. Recipe for constructing smooth wavelet frames using a single dilation 
matrix and an irregular grid. We can combine the results of this section with 
Theorem |3.1[ to obtain the following recipe to construct smooth, wavelet frames 
of M. d associated to a single dilation matrix A and an irregular grid. The wavelets 
obtained by this method can be constructed to have polynomial decay of any degree, 



as exemplified in Section 4.2 
Recipe: 

• Select a bounded set V C R d such that G V° and fi(dV) = 0. 

• Select a function h of class C r , r > such that \h\ ^ on Q = A T V \ V 
and Supp/i C Q E for some small e > 0. 

• Consider a set X = {xk}keK C K d , such that X is separated and p(X) < 
t^j, where 5 — diam(Q £ ). 

Then the following collection is a frame of L 2 (R d ): 

{\detA\ j / 2 ij{A j x-x k ) : keKJ e Z)}, 
where ip is the inverse Fourier transform of h. 



6. Reconstruction formulas 

We first no te th at although the set {| det Aj\ 1 ^ 2 ip(AjX — Xj^k) ■ j & J, k € K} 
in Theorem 3.1 is a wavelet frame for L 2 (R d ), it is not in general true that for a 
fixed j the set {tpj yXj k (x) = | det Aj\ 1 ' 2 ip(Ajx — Xj_k) '■ k £ K} is a frame, unless 
the wavelet frame is of Shannon-type, hence, not well-localized. Thus, for well- 
localized wavelets, it appears that the reconstruction of a function / € L 2 (R d ) 
from the wavelet coefficients {< /, "4>j, x - k >'■ j € J, k € K} cannot be obtained in 



a 

stable way by first reconstructing at each level j and then obtaining / by summing 
over all levels j. However what follows shows that this is still possible: 

For a matrix R, denote by Dr the dilation operator Dn(f) = \ det R^^foR. Using 
the notation of Theorem 13.11 we have that 



< /i^xj,* >=< f-l'il'i: < ,• . • ' /'.,/• n " ' ■ > ■ 
This, together with the assumptions that k XQ '■ k £ K} is a frame for JCq 
for each j, and that Supp/ij c Sj, permits to reconstruct fj = hjf using a dual 
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frame {4>j,k '■ k £ K} of {e Xj k XQ '■ k £ K}. Note that, since D B . is unitary then 
{D Bf <f>j,k ■ k £ K} is a dual frame of {D Bj {e. Xjtk Xo) ■ k £ K} for K Sj - 

Thus, it is always possible to reconstruct each fj in a stable way and then obtain 
/ by summing up over all levels j. One drawback is that the dual frame may not 
be well-localized, since each <f>j^ may be discontinuous at the boundary of Q. To 
treat this problem we simply note that 

hjf = ^2< f, hjD Bj e Xjik > D Bj <j>j,k- 
k 

Multiplying both sides by hj we obtain 



(6.1) 



N 2 



A,- 



where = h<pj,k- If we choose h to be in C r (M. d ), and if 4>j,k is m C r (Q) , 
then Oj.k will be in C r (M. d ) and therefore decays polynomially in space. Hence the 



partial sums of the series (6.1 1 will have good convergence properties. We can then 
sum equation (6.1) over j £ J to obtain fJ2j \hj\ 2 > an d then divide by \hj\ 2 
obtain /. 



to 
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